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We adopt methods that are well known in statistical physics to the problem of stochastic in- 
flation. The effective power spectrum for the classical, stochastic long-wavelength fluctuations is 
calculated for free scalar fields in a de Sitter background. For a smooth separation into long and 
short wavelengths, we identify an infra-red divergence of the effective power spectrum, which has its 
correspondence in statistical physics in the phenomenon of dimensional reduction. The inflationary 
dynamics pushes the affected scales exponentially fast to large superhorizon scales, and establishes 
scale-invariant behavior for smaller scales (for massless fields). In the limit of a sharp separation of 
wavelengths, the scale of the infra-red divergence is pushed to infinity. 
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The concept of stochastic inflation, introduced by 
Starobinsky several years ago [T], provides a framework 
to study the evolution of quantum fields in an inflation- 
ary universe |3] and has acquired considerable interest 
over the last years [HE]. The key idea lies in splitting the 
quantum fields into long- and short- wave length modes, 
and viewing the former as classical objects evolving in an 
environment provided by quantum fluctuations of shorter 
wavelengths. It constitutes an example of how the fun- 
damental properties of quantum fields can be modeled 
using methods of statistical mechanics. From this point 
of view, one focusses on the "relevant" degrees of free- 
dom, the long-wavelength modes, and regards the short- 
wavelength modes as "irrelevant" ones, generating a bath 
in which the former evolve. The natural length scale of 
this problem is the Hubble length, from which "short" 
and "long" acquire their physical meaning. 

This stochastic description of cosmological inflation 
has in fact two steps of complexity. In a simplified setup, 
the problem is reduced to describe a scalar test field (p 
in a fixed cosmological background. If ip is free, massive 
and minimally coupled, one obtains after splitting into 
long and short wavelengths, tp = ip L + tp s , an effective 
equation of motion of generalised Langevin-type, 

(□ + [i 2 )ip L (t 7 x) = h(t, x). (1) 

Here, p L is viewed as a classical entity, evolving stochas- 
tically in the presence of a (quantum) random force h, 
which is Gaussian distributed with zero mean. 

Early studies focussed on homogeneous fields — thus re- 
stricting attention to the time evolution of <p^. The study 
of inhomogeneous fields (see, e.g., [5]) is more involved, 
but also allows one to discuss spatial correlations. 

In the full problem of stochastic inflation, <p is the in- 
flaton field, rather than a test field in a fixed background. 
This requires not only a stochastic description of tp, but 
also of the geometry. Thus, one should analogously split 
the space-time metric and curvature into long- and short- 
wavelength parts. 

In this work, we study the scaling behavior and time 



evolution of the power spectrum of tp h in a fixed back- 
ground. We do so by means of replica field theory, which 
is well known in statistical physics [6 4 . As far as we know, 
replica fields have not been applied to cosmology before. 
They allow us to compute spatial correlations and their 
time evolution for the general test-field case. 

Below we illustrate this method by the example of 
a free, minimally coupled, ./V-componcnt, real field tp 
with mass /i. Let us further restrict our attention to 
a spatially-flat de Sitter universe in d space-time dimen- 
sions. Its scale factor is a(t) — e Ht with Hubble expan- 
sion rate H. For convenience we rescale to dimensionless 
variables tp/H^' 2 ^ 2 -> tp, fi/H -> fi, tH -> t, xH -> x, 
and k/H — > k, and use H = c = 1. The mode func- 
tion u(t, k) is defined via the decomposition of the field 
components, i = 1,...,N, 

<Pi(t,k) = ki(k)u(t,k) + h.c, (2) 

with k := \k\ and the annihilation and creation operators 
obey the commutation relations 

[&i(fc),at(p)] =(2n) d - 1 6 d - 1 (k-p)6 lJ , 
[ai(fe),aj(p)] = 0. 

The propagator is defined as 

G (t,t',k,k') ■= ^(n\<p(t,k)-<p(t',k')\n), (4) 

where the vacuum |f2) is defined by a,(k)\il) = at t = 
and a subscript "0" indicates a quantity that is calculated 
in the absence of any noise. 

An object of central interest in cosmology is the dimen- 
sionless power spectrum Vip(k). Its relation to some field 
propagator G(fc) with an infra-red behavior G(k) ~ k~ 71 , 
i.e., for k <C aH, is given by 

V v (k) := k^Gik) ~ k n -~\ (5) 

with the spectral index n v , connected to the critical ex- 
ponent rj via n v = d — rj. For d — 4 and /i = the power 



2 



spectrum of the free, noiseless theory is scale-invariant, 



1 for 770 = 3. Non-zero mass leads to 



, v = 4 - 3^/ 1 - -n 



2 _ 



= 1 



(6) 



Note that this result does not include any metric pertur- 
bation, which would cause the deviation from n v = 1 to 
be negative (see, e.g., [7]). 

Let us now introduce a split of the quantum field tp into 
short- and long- wavelength modes. For the free scalar 
field under consideration, this generates a linear noise 
term on the r.h.s. of the field equation, which becomes of 
type Q. On the level of the action, this corresponds to 
a linear random potential 



V(<£ L ) = h • <^ L 



(7) 



plus a term quadratic in h, which is ascribed to a Gaus- 
sian probability distribution p[{h}]. In statistical field 
theory, this is the so-called random-field case. 

For interacting scalar fields, the split of <p into long 
and short wavelengths generates higher powers of both 
<y9 L and the noise in the equation of motion. They may 
as well be collected into a probability distribution and a 
random potential, which might be written as 
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3=1 {<• - : 1 



(8) 



with the Taylor coefficients {h} being a set of random 
variables subject to p[{h}]. 

For any quantity O depending on {h}, the stochastic 
average, which shall be denoted by a bar, is then calcu- 
lated as 



0[{h}} := 2?[{h}]p[{h}]0[{h}]. 



(9) 



Note, that linearizing the equation of motion in the quan- 
tum modes, corresponds exactly to a Gaussian distribu- 
tion with vanishing mean. As has been shown recently, 
this effectively re-sums the leading- log contribution of the 
full quantum theory [?]. 

Objects of central interest are n-point correlation func- 
tions of the classical long- wavelength field ip L , i.e., 



<¥#(zi) ••••■¥#*(*»)>• (10) 
Since they can be derived from the generating functional 



Z[Ti] = J T>[<pj]exp{iS[<p h , j* L ]} 



(11) 



by taking derivatives of ln{Z[J^]} w.r.t. an external 
source J L , one in general needs to know ln{^[J L ]}. This 
might be difficult to calculate, because one has to aver- 
age a logarithm of a path integral over an exponential. 



To perform these stochastic averages we use the replica 
trick |U, 
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ln{Z} = lim — In \Z m \. 

m-+0 TO L ' 



(12) 



Thus, one just has to compute Z m for integer to, and, 
if the result is analytic in to, to take to — > at the 
end. This tantamounts to the introduction of to differ- 
ent copies (replicas) of the same system which are then 
coupled through the noise average. 

As already mentioned, the noise distribution p[{h}] is 
Gaussian for free scalar fields, with the first and second 
cumulants given by 



hi(x) = 0, 
hi(x)hj(y) = Sij<p(x,y), 



(13) 



with a known function cj)(x, y) depending on derivatives of 
the mode functions in ([2]) . In terms of the linear random 
potential ([7]), this may be written as 



V(ft«(aO) =0, 



(14) 



where the replica indices a,b = 1, . . . , to label the differ- 
ent copies arising by means of the replica trick. 

For interacting scalar fields, a Gaussian distribution 
is an approximation, but for many cases at least a rea- 
sonable staring point. Hence, one may specify a general 
random potential V as in QSJ) , in analogy to the case of 
free scalar fields, 



N 



V(<?S(x))V(0£(y))=4>(x,y)NR 



(15) 



where each of the random variables {h} is taken to be 
Gaussian-distributed with mean zero. Of course, the con- 
crete form of the function R, as well as the form of its ar- 
gument has to be determined from first principles. Note 
that for our example of a free scalar field, the Gaussian 
distribution is exact and that the function R is linear, as 



seen in ( 14 ) . In the terminology of statistical field theory, 
the space-time correlation cj>(x,y) is called short-range, if 
<t>{x,y) = p(t,t')d d ~ 1 (x — y), and long-range for all other 
cases (c.f. [9]). The latter is relevant for stochastic infla- 
tion. 

After averaging over the noise we obtain from Z m and 



equation (15) the replicated action 



S (m) = \ £ / / G L o\t, t\ k) #(t, k) ■ tf{t\ -k) 



m 

- J2 / ^y) NR 

L a,b=l Jx 'V 



N 



(16) 
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with the definitions J k := Jd^fe/^)^ 1 , j t := J dt, 
J x := J d^, and G L0 (fc) := J k , G L0 (fc, fc'). We dropped 
the source term oc j L , because it is irrelevant for the fol- 
lowing calculation. We see that the to replicas are cou- 



pled to each other through the noise correlation (15 1 



invariance for the spatial correlation, i.e., <f>(x,y) — 
4>{t, t', \x — y\). The function R is defined by R(( • ) ) := 
(R( • )) - In the limit N — > oo, or trivially for free scalar 
fields for arbitrary values of N, one has R = R. 

The physical interpretation of the saddle-point equa- 



In the following, we will use a variational method to tions (Il8| and ((20| is the following: The replica structure 



approximate action ( 16 I by the Gaussian variational ac- 
tion 



S, 



(m) 



- y 

2 ^ 

a, 6=1 



t.t' Jk 



G L - 1 a6 (t,t / ,fc)^„(t ) k)-0 b (t',-k) 



with the inverse propagator 



G 



l ab 



G 



lq Uab ^ab- 



(17) 



(18) 



Although Gaussian, the perhaps highly non-linear noise 
interaction R(... a f,) is accommodated by the non-diagonal 
replica structure a a b- 

Let us comment on the structure of ( 18 ). On the diag- 



onal (in replica space) we find the inverse of the noiseless 
propagator G L0 ~ 1 plus some mass correction, to be de- 
termined later. This correction alone would not only be 
trivial but also inconsistent, as we will see later. Hence, 
the off-diagonal part is filled by some, a priori unknown, 
replica structure a ab , which, in general, can be time de- 
pendent, and, if one includes long-range noise correla- 
tion, also momentum dependent, directly affecting the 
scaling behavior of the power spectrum. Thus, although 
this variational method only generates a self-energy con- 
tribution, its off-diagonal replica structure might have a 
viable influence on large-scale correlations. 

The Gaussian variational method becomes exact in the 
limit N — > oo and allows one to go beyond ordinary per- 
turbation theory. It is based on the following Fcynman- 
Jensen inequality |10j 



ln{Z}>ln{Z }+(4 m) -5(" l )) o , 



(19) 



where the subscript refers to the variational action (17) 



and we temporarily Wickrotate to Euclidean signature. 



Equation (19 1 can easily be proven by using the Jensen 
inequality exp{(. ..)} < (exp{. . .}) [11, which comes 
from the convexity of the exponential. The problem is 
to find the best G La 6, i-e., the best a a b, satisfying (19) 



by maximizing the r.h.s. of (19 1 



The result of the variation, again for Minkowski signa- 
ture, is 



(20) 

where cf>(t,x) := cf>(t,t,x), a ab (t,k) := a ab (t,t,k) and 
G Lo a(t, k) := G Lo j,(i, t, fc), the prime denotes a deriva- 
tive w.r.t. the argument, and we assume translational 



(j is a generalized self-energy (c.f. the discussion in sec- 
tion 3 of [6]). 



The replica structure ( 20 1 already contains our main 



result, namely the dimensional reduction at large scales, 
as will become clear in the remainder of this work. It 
also shows that the replica matrix is in general space- 
time dependent, which affects the scaling behavior of the 
two-point function (c.f. [12"]). As can easily be seen, in 
this (replica) Gaussian variational approximation, any in- 
teraction without random variables is diagonal in replica 
space. Therefore it only modifies the diagonal part of 
the self-energy. Arbitrary noise gives rise to a complex 
replica structure, however. 

For the free scalar field studied in this paper, one ob- 
serves that random-field noise (linear potential) neces- 
sarily yields a uniform replica matrix, i.e., a ab = cr for 
all a, b, because R' = const. This is the so-called replica 
symmetric case, which has intensively been studied [6]. 
It means that different replicas couple all in the same 
way among each other. To leading order in the number 
of replicas, to, we find that (G La f,) has the form 

(G Lafe )(i, k) = G L0 (t, k)l + a(t, k) G L0 (i, fc) 2 JJ. (21) 

where G L o(t, fc) := G L o(t, f, fe), and the m x m-matrix JJ 
has 1 in every entry. In the limit of vanishing correla- 
tion, i.e., a — > 0, the free theory is recovered. To derive 



(21 1, no assumption except the validity of the Gaussian 



variational method is made. 

The physical propagator G L of the long-wavelength 
field tp L is obtained from G La b via 

G L (i,fc)= lim -Tr[(G Lafc )(i,fc)], (22) 

m— >0 to 

which is simply the arithmetic mean of the trace of the 
replica matrix propagator (c.f. section 4 of (6]) and yields 



G L (t, fe) = G L0 (t, fc) + a{t, fe) G L0 (t, fe) 2 



(23) 



We now turn to the study of the infra-red behavior 
of the physical propagator and therefore the power spec- 
trum, with focus on the spectral index. Although irrele- 
vant for stochastic inflation, we first assume short-range 



correlation for pedagogical reasons. Equation (20 1 then 
implies a to be independent of momentum. With the 
infra-red behavior G L0 (i, fe) ~ fc _??0 and G h (t,k) <~ k~ v 



(770,77 > 0), we obtain with the help of (211 and (23 1 the 
relation 



77 = 2?7o 



(24) 
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— drastically different from to the noiseless case, where 
7/ = 7/o- This is a variant of the phenomenon of dimen- 
sional reduction |T3l [14]. It is best understood in x- 
space, where the two-point function goes in the infra-red 
like \x\ ^ 1 - )+T ' 5 and the change in the exponent 77 can 
be absorbed into the reduced dimension d! := d — 7/ . 
We have to stress that, in the random-field case and for 
short-range correlation, the dimensional reduction theo- 
rem can rigorously be proven to hold in great generality, 
i.e., to all orders in perturbation theory and for arbitrary 
non-random potentials (see especially [H] for a super- 
symmetric version of the proof). 

Let us now consider the case of long-range correlation 
of the form <p(x) ~ \x\ ^ l ' +p with p < d — 1. This de- 
scribes properly the infra-red limit of the physical model 
discussed below. In momentum space, the above choice 
implies <p(k) ~ k~~ p and hence a a b(k) ~ k~ p by virtue of 
(20 1. For p > —770, the infra-red behavior of the power 



spectrum deviates from the noiseless result, and we find 
77 = 2770 + p. This result is consistent with previous stud- 
ies in flat space with 77 = 2 [T2"] . 

For exponential inflation (d = 4 and p = 0) we find a 
modification if p > —3, i.e., the spatial noise correlator 
decreases at most like \x\ 6 . In general, the spectral 
index changes to 



d-2r] Q - p, 



(25) 



wherefore the choice p = yields the short-range result. 

The consequence is a dramatic change of the super- 
horizon power spectrum of the classical long- wavelength 
modes as compared to the case without noise. In partic- 
ular, scale invariance does not hold on these scales. 

Let us now return to our physical model of stochastic 
inflation. The split of the field <p into a long- and short- 
wavelength part, tp — ip^ + ips, together with the free field 
equation, (□ + p 2 )ip — 0, implies 



a(t,k) = 



W, 



a(t) 



— e I u(t, k) 



(26) 



where is the (spatially Fourier transformed) covariant 
Laplacian, u(t, k) is the mode function from equation ([2]), 
and W K is a smooth high-pass filter, cutting out the low 
frequencies below e. The parameter k controls the width 
of the cut. In the limit n — > 0, W K approaches a step 
function. Here we choose W K ( • ) = l/ir arctan( • /k) +1/2 
and take < e <C 1 in order to cut far below the Hubble 
rate H(= 1), and k <C e to have a narrow transition 
region between quantum and classical modes. We do not 
impose any restriction on p except that we demand the 
radicant in Q to be positive, i.e., p 2 < 9/4. 

It is natural to define a transition scale fc* at which 



the two terms on the r.h.s. of equation (23) balance each 
other. It separates two regions such that for k 3> fc* 
the behavior is noiseless and for k <ti fc*, dimensional 
reduction holds. 



In the infra-red, the model of equation (26) is of long- 
range type. We find in four dimensions 



and thus with ( 25 ) for k <C k 



6-9W1 
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-3 + 2/i 2 + 0(p 4 



(27) 



(28) 



while for k 3> k* the noiseless spectral index Q is recov- 
ered. 

The zeros of the relative correction 



G L (*,fc)-G L0 (M) 
G L o(i) k) 



ff(i,A)G L0 (t,fc) (29) 



define the transition scale k*(t). Its late-time behavior 
can be calculated analytically, 



2-\/9-4^ 2 -2 



8-2^9 1/.- 



(e 2 + K 2 )\ 
For p = we find the asymptotic form 



) 



^tt\/ e 2 + k 2 



(30) 



(31) 



Thus, for £ / 0, i 4 goes to zero in the (step-function) 
limit k — > 0, i.e., dimensional reduction disappears. We 
should point out that the two limits k — * and e — * 
do not commute and that the case k ^ and e — * is 
unphysical. 
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Figure 1: Comoving transition scale fc* as a function of 
cosmic time t (both in units of H) for mass p/H = 
0.1, 0.2, 0.3, 0.4, 0.5 and 0.6 (dashed lines, top to bottom). 
Dashed curves are numerical results, colored solid lines are 
analytic approximations, and enveloping black lines are | a(t) 
and the asymptotic form (311, respectively. A smoothing 
K — 10 -3 and short-wavelength cut e = 10 -2 are chosen. 

Figure [T] shows the time behavior of the comoving 
scale fc* for different values of the mass p. The solid 



5 



rays represent the analytic approximation ( 30 ) , while the 



dashed curves are obtained numerically from the roots of 



( 29 1 . Well below this borderline the two-point function of 



the classical stochastic field obeys dimensional reduction, 
while well above ordinary scaling holds. One sees that 
fi — gives the largest k*(t). After an initial transient 
phenomenon, whose duration depends on the specific 
choice of e and n, the comoving transition scale decays ex- 
ponentially fast. Hence, the dimensional- reduction con- 
tribution is pushed to larger and larger scales as time 
increases. This therefore guarantees that quantum noise 
induces only a minor change of the spectral index on sub- 
horizon as well as on moderate super-horizon scales. 

For concreteness, let us consider a mode with comov- 
ing k = 0.05iJ. At time t = it is within the region of 
ordinary scaling, suffering at most slightly from dimen- 
sional reduction. This mode enters then, after roughly 
two e-foldings, the region of broken scale invariance, but 
leaves it at the latest (for /i = 0) after seven e-foldings 
and stays eternally in the scale-invariant regime, which 
itself grows exponentially fast. 

One may connect the replica structure a to a non- 
linearity parameter g NL , which shall now be defined via 



<Pi(t, k) = cpf{t, k) - g NL (t, k)(<pf(t, k))' 



(32) 



where ip G (t,k) is a free Gaussian field. On the level of 
propagators, this translates to 



G L (t, k) = G L0 (t, k) + 3 9nl (*, k) 2 G L0 (t, kf 



and hence 



a(t,k) = 3g NL (£, kf 



(33) 



(34) 



can be directly read off, using equation (23 1. g NL mea- 



sures the influence of the quantum fluctuations, picked 
up by a smooth filter function. Formally, it resembles 
an effective non-Gaussianity parameter |15j for the long- 
wavelength modes. However, this association is mislead- 
ing since the theory we work with is Gaussian (but with 
non-trivial replica structure). 

Figure [2] shows the dependence of g NL on the comov- 
ing momentum k for various values of [i for fixed time 
t = 10H, using equation (34 1. Firstly, one sees that 
increasing /i lifts the curve upwards, and secondly, one 
observers a divergence in the infra-red — displaying the 
effect of dimensional reduction. For k 3> k* one obtains 
a scale-invariant spectrum. 

To summarize, we presented novel methods to study 
the power-spectrum of classical, stochastic fields in 
stochastic inflation. We demonstrate that replica field 
theory allows us to study the spatial behavior of non- 
coincident long- wavelength correlation functions. Di- 
mensional reduction changes the spectral index on super- 
horizon scales, in the sense that it heavily amplifies the 
power spectrum of the classical modes in the infra-red. 
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Figure 2: Non-linearity parameter p NL as a function of co- 
moving momentum k (in units of H) for mass /J./H = 
0.1 (lowermost), 0.2, 0.3, 0.4, 0.5 and 0.6 (uppermost) with 
Ht — 10, k — 10" 3 , and t = 10~ 2 . The dotdashed black 
line represents the value of fc* corresponding to the respective 



This effect has been calculated by a variational method, 
which allows us to go beyond ordinary perturbation the- 
ory and suggests it might be entirely of non-perturbative 
nature. However, in the limit of vanishing k, i.e., of a 
sharp cut between long- and short- wavelength modes, the 
singularity disappears. 

Those huge differences in the energy density on large 
super-horizon scales signal a breakdown of ordinary per- 
turbation theory, because in that case one cannot speak 
about the spectrum of fluctuations in the usual, pertur- 
bative sense. It further displays the failure of the test- 
field assumption, since in the situation at hand it is no 
longer valid to neglect the back-reaction of the field on 
the geometry. However, the time-evolution of the long- 
wavelength field pushes the dimcnsionally reduced region 
exponentially fast to unobservable scales. This provides 
further support for the self-consistency of the idea of in- 
flation, since regions of broken scale-invariance, with ex- 
traordinarily large fluctuations disappear faster than any 
causal patch of the universe expands. 

It will be interesting to discuss other background 
space-times, as well as to include self-interactions. The 
present formalism might also help to tackle the full prob- 
lem of stochastic inflation, where the geometry is also 
random. If we could assume that the dimensional reduc- 
tion of our test field also applies to the inflaton, we would 
find further support of the idea of eternal inflation [5] . 
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